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Abstract 

The generalized nonlinear Klien-Gordon equation plays an important role in quantum mechanics. In this paper, a new three- 
time level implicit approach based on cubic trigonometric B-spline is presented for the approximate solution of this 
equation with Dirichlet boundary conditions. The usual finite difference approach is used to discretize the time derivative 
while cubic trigonometric B-spline is applied as an interpolating function in the space dimension. Several examples are 
discussed to exhibit the feasibility and capability of the approach. The absolute errors and Lo^error norms are also 
computed at different times to assess the performance of the proposed approach and the results were found to be in good 
agreement with known solutions and with existing schemes in literature. 
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coUisionless plasma, the recurrence of initial states and examina- 
tion of the nonlinear wave equations was in [1]. 

Several methods, in addition to several finite difference schemes, 
have been developed to solve the nonlinear KG equation. Jiminez 
and Vazquez [2] introduced four numerical schemes for solving 
the nonlinear KG equation. Ming and Guo utilized a Fourier 
collocation method for solving the nonlinear KG equation [3]. 
The KG equation was approximated using decomposition scheme 
by Deeba and Khuri [4] and using the Legendre spectral method 
by Guo et al [5]. Wong et. al. solved an initial value problem 
involving the nonlinear KG equation using fuUy implicit and 
discrete energy conserving finite difference scheme [6]. Wazwaz 
introduced the tanh and sine-cosine method to obtain compact 
and noncompact solutions for the nonlinear KG equation [7]. 
Sirendaoreji solved the nonlinear KG equation using the auxiliary 
equation method to construct new exact traveling wave solutions 
with quadratic and cubic nonlinerity [8]. Yucel solved the 
nonlinear KG equation using homotopy analysis method [9] 
and Chowdhury and Hashim solved the equation using homo- 
topy-pertubation method [10]. 

B-spline functions can be used to solve numerically linear and 
non-linear differential equations. Caglar et. al. [1 1] has introduced 
a cubic B-spline interpolation method to solve the two-point 
boundary value problem. Hamid et al. [12] has introduced an 
alternative cubic trigonometric B-spline interpolation method to 
solve the same problem. Dehghan and Shokri [13] have obtained 
a numerical solution of the nonlinear KG equation using Thin 
Plate Splines radial basis functions. IGiuri and Sayfy [14] have 
solved the generalized nonlinear KG equation using a finite 



Introduction 

The generalized nonlinear Klien-Gordon (KG) equation arises 
in various problems in science and engineering. This paper focuses 
on the analysis and numerical solution of the generalized 
nonlinear KG equation, which is given in the following form [14]: 

^u{x,t) + cc^u{x,t)+liu{x,t) + Giu{x,t))=f{x,t), 
a<x<b, 0<t< T 
subject to initial conditions 

{u{x,0) = wi (x), 
8 a<x<b (2) 

— u{xfi)=0)2{x), 

and with Dirichlet boundary conditions 

' u{a,t) =&'. (t), , , 

^ ' ' 0<t<T 3 

where u = u{x,t) denotes the wave displacement at position and 
time (x,f), a and ji are real constants, G{u) is a nonlinear function 
of M and f{x,t),03\(x),032{x),^\{t) and ^ji^) '^'"e known 
functions. 

In particular, the KG equation is important in mathematical 
physics especially in quantum mechanics and it is well known as a 
soliton equation. A study on the interaction of soliton in 
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Figure 1. (a) Analytical solution, u(xj) of Problem 1 (b) Approximate solution, u{x,r) of Problem 1 with ii = 40 and A^^O.l. 
doi:10.1371/journal.pone.0095774.g001 



element coUocation approach based on third degree B-spline 
polynomials. 

In this work, a new three-time level implicit approach which 
combines a finite difierence approach and cubic trigonometric B- 
spline collocation method (CTBCM) is proposed to solve 
generalized nonlinear KG equation. The finite difFerence 
approach is proposed to discretize time derivative and cubic 
trigonometric collocation method is appKed to interpolate the 
solutions at time/. Two numerical experiments are carried out to 
calculate the numerical solutions, absolute errors, Loo error norms 
and order of convergence for each problem in order to show the 
accuracy of method. 

Temporal Discretization 

Consider a uniform mesh Q with grid points to discretize 

the grid region A=[a,h] x [0,T] with Xj=a+jh, j = Q,\,l,...,n and 



tk = kht, k = Q,\,2,l>,...,N , r = /i and Ar denote mesh space 
size and time step size respectively. The time derivative is 
approximated using the central finite difference formula 



M? ' '-lu'f + u''- 



(4) 



Using the approximation of equation (4), equation (1) becomes 



„A'+1 



(AO 



+ a^^Pi^=f[xj4i,)-G[u[xjJk}) (5) 



Analytical and Approximate Solution 
0.6 1 ■ . . 1 1 




Figure 2. (a) Analytical and approximate solution at different time (b) Error of Problem 1 with /i = 40 and A^^O.l. 
doi:1 0.1 371/journal.pone.0095774.g002 
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Table 1. Numerical solution of Problem 1 at ? = 5 





xjn 


0 


0.2 


0.4 


0.6 


0.8 


1.0 


10 


0 


0.0567 


0.1134 


0.1702 


0.2269 


0.2837 


20 


0 


0.0567 


0.1135 


0.1702 


0.2269 


0.2837 


40 


0 


0.0567 


0.1135 


0.1702 


0.2269 


0.2837 


80 


0 


0.0567 


0.1135 


0.1702 


0.2269 


0.2837 


u{x,t) 


0 


0.0567 


0.1135 


0.1702 


0.2269 


0.2837 



doi:1 0.1 371 /journal.pone.0095774.t001 



Using the 6— weighted technique, the space derivatives of 
equation (5) becomes 



,k+\ , 



-2u'' + u' 



{At) 



^ + (i-%;+%|+'=y;*^-G(«;) (6) 



and the subscripts k and k+l are successive time levels 
After simplification, equation (6) leads to 



P{u)j+\ o<e<i 



{i-eKA( 



(7) 



Trigonometric B-Spline Collocation method 

In this section, CTBCM is used to solve Klien-Gordon 
equation. Cubic trigonometric B-sphne are used to approximate 
the space derivatives. To construct the numerical solution, nodal 
points (Xj,ti^) defined in the region [a,b] x [0,2"] where 
a = X(i<xi< ■■ ■ <x„ = h, Xj^i—Xj = h. 

The approximate solution u{xj) to the exact solution u{x,t) is 
defined as [18]: 



p{xj){p(xj)q(xj+2) + q(xj+3)p{xj+i))+q(xj+i)p^{xj+i), xe[xj+uXj+2] (9) 
'!ixj+i){p{Xj+i)qixj+i) + q{xj+i)p(xj+2))+p{xi)q^(xj+,), xe[Xj+2,Xj+,] 
, xe[xj+j,Xj+i\ 



where 



q(xj)= sin(-^y^), 

h y-i 

K= sin(-)sin(/2)sin(y). 



Due to local support properties of B-spline basis function, there 
are only three non-zero basis functions 74j_3(x/), T^j^iiXj) and 
T^j-llXj) are included over subinterval + Thus, the 

approximation u'- and its derivatives with respect to x can be 
simplified as: 



u(xj)= Ci(mj(x) 



(8) 



where Cj{t) are time dependent unknowns to be determined and 
T4j(x) are cubic trigonometric B-spline basis function given as: 



where 



= ;?i C/_3(0 + ihCj-iit) + niCj-i(t) 
{ux)j = -ri3 Cj-3(t) + % Cj- 1 (0 

{Uxx)j = tl4 Cj-l(t) + C,_2(0 + (74^,-1 (0 



))i = csc(/i)csc( y I sm ( - 1, 



(10) 



Table 2. Absolute error of Problem 1 at ? = 5. 





x/n 


0 


0.2 


0.4 


0.6 


0.8 


1.0 


10 


0 


2.473 X I0-= 


3.810X 10-5 


3.559 X 10-5 


2.152X 10- = 


0 


20 


0 


6.486 X IQ-'' 


1.496 X IQ-' 


8.892 X 10-" 


5.485 X 10-« 


0 


40 


0 


1.709 X lO-"* 


2.490x10-' 


2.337 X lO-" 


1.442 X 10-' 


0 


80 


0 


5.119X IQ-' 


7.460 X 10-' 


7.003 x 10-' 


4.320 X 10-' 


0 



doi:1 0.1 371 /journal.pone.0095774.t002 
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Table 3. Comparison of L^o errors norms with Khuri & Sayfy [14] using h = OA, A7 = 0.005 and Dehghan and Shokri [13] using 
/! = 0.02, A? = 0.0001 





t 


1 


2 


3 


4 


5 


Dehghan& Shokri [13] 


1.254 X 10-5 




1.555 X 10-5 




3.379 X 10-5 


Khuri& Sayfy [14] 


2.838 x 10-'' 


3.299 X 10-* 


7.005 X 10-5 


3.018 X 10-* 


3.249 X 10-* 


Present method 


4.552 X 10-5 


4.358 X 10-5 


1.359X 10-5 


5.353 X 10-5 


3.868 x 10-5 



doi:l 0.1 371 /journal.pone.0a95774.t0a3 



12-- 



>h-- 



'74 = 



l+2cos(/i) 

3 /3h 

- tCSC — - 

4 V2 
3 + 9cos(/i) 



4cos(/!/2)-4cos(5V2)' 
3cot2(/!/2) 



'75 = 



2 + 4cos(/z) ' 



The solution to problem (1) is obtained by substituting equation 
(10) into equation (7). Initially, time dependent unknowns C" are 
calculated and shown in section 3.1. Next, the following initial 
condition is substituted into the last term of equation (7) for 
computing C' 



= Uj —2Ata)2{x) 



(10a) 



Subsequendy, the time dependent unknowns, C'' for A: > 1 
should be generated. After simplification, equation (7) leads to 



(11) 



where. 



P5 = 1 + 0yi, 
P3 = 2-(l-%i, 



with 



P4 = (l-%2, 

d]=D{xj4k)=(^tf(Ff-G{j] 
y2 = ci{Atf. 

The system obtained on simplifying (11) consists of w + 3 
unknowns C'' = (C*^3,C^2;C'* j, ■ ■ ■ ,C*_[) in «+l linear equa- 
tions at the time level f = r/,-+i. In order to obtain a unique 
solution, the equation (8) is apphed to the boundary conditions 
given in equation (3) for two additional linear equations. 

u{a,tk+i) = riiCj-3(t) + ri2Cj-2(t) + ihCj-i(t) = (tii{tk+i) (12) 
u{b,tk+i) = r]iCj-^(t) + ri2Cj-2(t) + riiCj-i(t) = <l>2{tk+i) (13) 

From equations (1 1), (12) and (13), the system can be written as 



MC''+'=NC*-PC*"'- 



(14) 



where. 



Table 4. The maximum errors norms 


and order of convergence, 


p for Problem 1. 






' = 2 




t = 5 










P 


10 4.358 x 10-5 




3.868 x 10-5 




20 1.094x10-5 


1 .9942 


9.643 x 10-0 


2.0039 


40 2.886x 10-' 


1.9927 


2.532 X 10-' 


1 .9290 


80 6.648x 10-' 


1.7384 


7.592 X 10-' 


1.7379 


doi:l 0.1 371 /journal.pone.0095774.t004 
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Table 5. Numerical solution of Problem 1 at ? = 5 





xjn 


0 


0.2 


0.4 


0.6 


0.8 


1.0 


25 


0 


0.0567 


0.1134 


0.1702 


0.2269 


0.2837 


50 


0 


0.0567 


0.1135 


0.1702 


0.2269 


0.2837 


100 


0 


0.0567 


0.1135 


0.1702 


0.2269 


0.2837 


200 


0 


0.0567 


0.1135 


0.1702 


0.2269 


0.2837 


u{x,t) 


0 


0.0567 


0.1135 


0.1702 


0.2269 


0.2837 



doi:l 0.1 371 /journal.pone.0095774.t005 



M 



0 
0 



N = 



1 




ni 


0 


• 0 


0 


0 0 




M2 Mt, 


0 






0 




Ml M2 


Ml, 






0 










M\ 


M2 


M3 0 










0 


Ml 


M2 M3 




0 


0 


0 


■ 0 


n\ 


n2 n\ 


/o 


0 0 


0 


■■ 0 


0 


0 0 


\ 




N2 N3 


0 






0 




0 


Ni N2 


iV3 






0 




0 






TV, 


N2 


7V3 0 




0 






0 


Ni 


N2 N3 




0 0 


0 


■■ 0 


0 


0 0 


/ 




/o 


0 0 


0 


■■ 0 


0 


0 0 \ 






'7i 




0 • 






0 






0 


n\ n2 


'7i 






0 






0 






m 




rji 0 






0 






0 




'72 '71 






U 


0 0 


0 


■ ■ 0 


0 


0 0/ 





With 



{Atf[F^-G{u'i)] 



(ao'K-gK) 



Mi=Pini+P2{-%), 

M2=Pslh, 
M3=P5tll+P2{th), 

Ni=Pin\-p^{-m), 

^2=P3'72' 
A^3=P3'7l-P4('73)- 



Half explicit and half implicit scheme is produced by choosing d 
to be 0.5. This scheme is known as Crank-Nicolson scheme. 
System (14) becomes a tri-diagonal matrix system of dimension 
(m + 3) X (« + 3) that can be solved using the Thomas Algorithm 
[17]. 



Table 6. Absolute errors of Problem 1 at ? = 5. 





xjn 


0 


0.2 


0.4 


0.6 


0.8 


1.0 


25 


0 


4.404 X 10-'' 


6.172X I0-' 


5.745 x 10-'' 


3.524 X 10-<* 


0 


50 


0 


1.101 X 10-" 


1.550X 10-' 


1.439 X 10-'' 


8.827 X 10-' 


0 


100 


0 


2.783 X 10-' 


3.926x10-' 


3.644x10-' 


2.236 X 10-' 


0 


200 


0 


7.320 X 10-* 


1.033 X 10-' 


9.587 x 10-* 


5.882 X 10-* 


0 



doi:l 0.1 371 /journal.pone.0095774.t006 
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Table 7. Comparison of Loo errors norms with Khuri & Sayfy [14] using h = 0.04, A? = 0.001 and Dehghan and Shokri [13] using 
/! = 0.02, A? = 0.0001 



t 


1 


2 


3 


4 


5 


Dehghan& Shokri [13] 


1.254 X 10-5 




1.555X 10-5 




3.379 X 10-5 


Khuri& Sayfy [14] 


4.599 X 10-5 


8.053 x 10-5 


1.276X 10-5 


7.292 X 10-5 


5.128x10-5 


Present method 


7.316X 10-' 


6.986 X lO-*" 


2.089 X 10-"* 


8.596 X 10-0 


6.245 X 10-0 



doi:l 0.1 371 /journal.pone.0095774.t007 



Initial vector C" 

The initial vectors C^can be obtained from the initial condition 
as well as boundary values of the derivatives of the initial condition 
[11,15]: 



u]=rnCl,+n2Cf_^+,nC^_^ = wi{xj\ y=o,i,2,...,« (i5) 

_ («0),= -,,3CO_3 + ,73C;Li =ir;(x,), j = n 

This yields a (« + 3) x (n + 3)matrix system: 
AC'' = B 



/-% 0 ,,3 0 

0 j;, t)2 (J, 

0 
0 

\0 0 0 0 



0 0 0 0 \ 
0 
0 



0 n 

0 0 nJ 



Ml (a-o) 

coi (x„) 



(16) 



The solution can be obtained by using the Thomas Algorithm 
[17]. 

Stability analysis using von Neumann method 

The von Neumann analysis of stability considers the growth of 
error in a single Fourier mode [19-20] 



Cf = S''exp{iiijh) 



(17) 



where / = v — 1 and f/ is the mode number. It is known that this 
method can be used to analyze the stability of linear scheme. All 
the nonlinear term in (11) are assumed to be zero [19]. Thus, 
equation (11) becomes 



-(l-%2("xx)--«f 



(18) 



where >'[ = {At) (i and 72 = (^0 Substituting the approximate 
solution (8) in equation (18) leads to 



where 



-'hCj_2 -liCj^i 
Pi =1i+S{Pi), 

P2=ni + 6{P2), 

P3=2;?, + (l-e)(p,), 
P4=2;j2 + (l-e)(p2), 

Pi='y4yi-^'?i72: 
P2='?5yi+'?2y2- 



(19) 



Table 8. The maximum L^o error 


norms and order of convergence, p 


for Problem 1. 






' = 2 




t=i 






p 




/' 


25 6.9860x 10-' 




6.2453 x 10-' 




50 1.7540x 10-' 


1.9938 


1.5704x10-' 


1.9916 


100 4.4432x 10-' 


1.9810 


3.9775x 10-' 


1.9812 


200 1.1690x 10-' 


1 .9263 


1.0466 X 10-' 


1 .9262 


doi:l 0.1 371 /journal.pone.0095774.ta08 
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(a) 



(b) 




0 0 




Figure 3. (a) Analytical solution, u{x,t) of Problem 2 (b) Approximate solution u(x,t} of Problem 2 with n = 5, Ar = 0.01. 
doi:1 0.1 371 /journal.pone.0095774.g003 



In order to obtain the amplification factor \5\, the trial solution 
(17) is substituted in (19) and after some simplification, we obtain, 



5'^A + SB+C = 0 



(20) 



where A=pi{2cosiih) +p2, B =pT,{2 cos tih)+p4 and 
C = t^i(2cosiih) + ri2- Since A,B,C>0 and 0<9<1, the ampli- 
fication factor of this scheme is 



Numerical Results and Discussions 

In this section, the CTBCM is applied on two numerical 
problems. In order to measure the accuracy of the method, 
absolute errors and Loci error norms are calculated using the 
following formulas [16] 



Absoluteerror = — Uj\ 



(22) 



\d\- 



M 



M + ON 



<1 



(21) 



where M = )7j (2cos f;/;) + );2 and = P[ (2cos (;/;) +P2- Hence, 
this scheme is unconditionally stable. 



max M, — w, 



(23) 



where and M,- are analytical solution and approximate solution of 
proposed problem (1), respectively. The numerical order of 
convergence, p is obtained by using following formula [14] 
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Table 9. Numerical solution of Problem 2 at ? = 1 





xjn 


0 


In 2/5 


21n2/5 


31n2/5 


41n2/5 


In 2 


5 


3.1416 


3.1222 


3.0654 


2.9733 


2.8495 


2.6990 


10 


3.1416 


3.1223 


3.0655 


2.9734 


2.8496 


2.6990 


20 


3.1416 


3.1223 


3.0655 


2.9734 


2.8496 


2.6990 


40 


3.1416 


3.1223 


3.0655 


2.9734 


2.8496 


2.6990 


u{x,t) 


3.1416 


3.1224 


3.0657 


2.9736 


2.8497 


2.6990 



doi:1 0.1 371 /journal.pone.0095774.t009 



_ Log(La, («)) - Log(Loo (In)) 

where [n) and Lqo (2w) are the Z-oo at number of partition n 
and 2n respectively. 

Problem. 1 

Consider the following nonlinear Klien-Gordon equation 
[13,14] 

u„ — Mv.Y -\-u^ = —X COS t + cos^ ; (25) 
subject to the following boundary and initial conditions 

u{xfi)=x, u,{x,Q)=Q; 0<.\:<1 



problem is calculated by the use of the formula given in (24) and is 
tabulated in Table4 and TableS. An examination of these tables 
indicates the method has a nearly second order of convergence. 

Case.l. Numerical solutions, absolute errors, error norms 
and order of convergence using time step size A/ = 0.005 

Case.2. Numerical solutions, absolute errors, L^i. error norms 
and order of convergence using time step size A? = 0.001 

Problem.2 

The following nonlinear Klien-Gordon equation which is also 
known as Sine-Gordon equation is considered [14] 

iiti — Uxx + sin M = 0 (26) 

It is subjected to initial conditions and boundary conditions as 
«(x,0)=0, M,(x,0)=4sechA-; 0<a-< ln2 



u(Q,t) = Q, u{l,t)= cost; 0<t<5 

The analytical solution of problem (25) is known to be 
u{x,t)=xcost and graphically shown in Fig. 1 (a). This problem 
is tested using different values oth and At to show the capability of 
the present method. The final time is taken to be T = 5.0. Fig. 1 (b) 
and Fig. 2 (a) show the approximate solution and Fig. 2 (b) shows 
the error of this problem with n = 40 and A/ = 0.1. Two cases of 
this problem are discussed where Case 1 and Case 2 consider 
A/ = 0.001 and A? = 0.005,respectively. Numerical solutions, 
absolute errors and order of convergence of each case are 
tabulated in Tablesl-4 and Tables5-8. The LooCrror norms are 
compared to Dehghan and Shokri [13] and Khuri and Sayfy [14] 
in Table3 and Table?. The comparison indicates that the present 
method is more accurate. The order of convergence of the present 



M(0,?)=4tan"' u{\n2,t)=4tw^ 0<r<l 

The analytical solution of this problem is 
M(x,r)=4tan^'(?sechx). Fig.3 (a) depicts a graph of this 
analytical solution. The final time is taken as7' = 5.0. The 
approximate solutions are calculated at time step size 
A? = 0.01with different mesh space size, h. Numerical solutions 
are recorded in Table9 and graphical solutions are plotted in Fig.3 
(b) and Fig.4 (a). Absolute errors are calculated and shown in 
Table 10 while the 3D error plot is depicted in Fig.4 (b). Table 11 
shows the comparison of Lao error norms between the present 
method with Khuri and Sayfy [14] method. This comparison 
shows that the present method gives better results. 



Table 10. Absolute error of Problem 2 at ? = 1 



xjn 


0 


In 2/5 


2 In 2/5 


3 In 2/5 


4 In 2/5 


In 2 




5 


0 


2.686 X 10-" 


3.534 X IQ-* 


3.072 x 10-* 


1.942 X 10-* 


0 


10 


0 


1.744 X 10-* 


2.525 x 10-* 


1.994 x 10-* 


1.080 X 10-* 


0 




20 


0 


1.492 X 10-* 


2.266 X 10-* 


1.737X 10-* 


8.401 X 10-' 


0 


40 


0 


1.435 x 10-* 


2.193x10-* 


1.678 X 10-* 


7.786 X 10-' 


0 



doi:l 0.1 371 /iournal.pone.0095774.t01 0 
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A New CTBCM to Numerical Solution of Generalized Nonlinear KG Equation 
Table 11. Comparison of Loo errors nornns with Khuri and Sayfy [14] using // = (ln2)/5, Ar^O.Ol 



t 0.01 0.02 0.1 0 0.50 1 .00 

Khuri& Sayfy [14] 1.3 x 10"^ 6.6 x 10"^ 3.8 x 10"'' 4.8 x 10"^ 2.2 x 10"^ 

Present method 5.566x lO"'' 1.618x 10"^ 9.851 x 10"^ 1.323x 10"-^ 3.534x lO"'^ 



doi:1 0.1 371 /journal.pone.0095774.t01 1 

Conclusions 

In this work, Klien-Gordon equation has been successfully 
solved using CTBCM incorporating a finite difference scheme. 
Specifically, the central difference approach is used to discretize 
the time derivatives and cubic trigonometric B-spline is used to 
interpolate the solutions at displacement x. Well-known two test 
problem were solved using the proposed method and the solution 
obtained were in good agreement with the known solution. 
Accurate solutions at intermediate points can be easily obtained. 
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